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Abstract 

Context. X ray clusters are conventionally divided into two classes: "cool core" (CC) and "non cool core" (NCC) objects, on the basis 
of the observational properties of their central regions. Recent results have shown that the cluster population is bimodal (Cavagnolo 
etal. 2009). 

Aims. We want to understand whether the observed distribution of clusters is due to a primordial division into two distinct classes 
rather than to differences in how these systems evolve across cosmic time. 

Methods. We systematically search the ICM of NCC clusters in a subsample of the B55 flux limited sample of clusters for regions 
which have some characteristics typical of cool cores, namely low entropy gas and high metal abundance 

Results. We find that most NCC clusters in our sample host regions reminiscent of CC, i. e. characterized by relative low entropy gas 
(albeit not as low as in CC systems) and a metal abundance excess. We have dubbed these structures "cool core remnants", since we 
interpret them as what remains of a cool core after a heating event (AGN giant outbursts in a few cases and more commonly mergers). 
We infer that most NCC clusters have undergone a cool core phase during their life. The fact that most cool core remnants are found 
in dynamically active objects provides strong support to scenarios where cluster core properties are not fixed "ab initio" but evolve 
across cosmic time. 

Key words. Galaxies: clusters: general - X-rays: galaxies: clusters 



1. Introduction 

Galaxy clusters are often divided by X-ray astronomers into two 
classes: "cool core"(CC) and "non-cool core" (NCC) clusters. 
The former are characterized by a set of properties: a promi- 
nent surface brightness peak, usually roughly coincident with 
the center of the large scale X-ray isophotes and with the posi- 
tion of the brightest central galaxy (BCG), associated to a de- 
crease of the temperature profile in the inner regions and to 
a positive gradient of the metal abundance profile. In the cen- 
tral regions of these clusters, the cooling time is significantly 
smaller than the Hubble time but the radiative losses are bal- 
anced by some form of heating, that in the currently prevailing 
scenario is attributed to the central AGN, which prevents the 
gas from cooling indefinitely and flowing inward. NCC clus- 
ters are characterized by the lack of these observational fea- 
tures. Several indicators based on these observational charac- 
teristics have been proposed to classify clusters into these cat - 
egories: e.g. temperature drop (Sanderson et al.ll2006l |2009b). 
cooling time (Bau er etldl |2005|). a composition of these two 
criteria dDunn et alJ 12005b iDunn & Fab ian 2008 ;), the slope of 
the gas density pro file at a given radius ([yikhlinin et al.l f2007) 
and core entropy dCavagnolo et al.l 2009). In a recent paper 
dLeccardi et alj|2009l hereafter Paper I), we suggested a robust 
indicator to classify clusters (see also Sect. [3j, based on pseudo- 
entropy gradients. The classification based on this indicator im- 
proves the traditional classification scheme based on the temper- 
ature drop and is essentially equivalent to classifications based 
on the cooling time. Increasing attention has been devoted to the 
statistics of CC, both in the local Universe and at higher red- 
shifts. While the precise results depend strongly on the indica- 
tor used to classify clusters, the fraction of CC is considered to 
be about 50% of the clusters population. There are some indi- 



cation that the cluster population is bimodal (Cavagnolo et al. 
2009) but there are also some intermediate objects which are not 
easily classified (Paper I). 

One of the open questions in the study of galaxy clusters, con- 
cerns the origin of this distribution. The original model which 
prevailed for a long time assumed that the CC state was a sort 
of "natural state" for the clusters and the observational features 
were explained within the context of the old "cooling flow" 
model: radiation losses cause the gas in the centers of these clus- 
ters to cool and to flow inward. Clusters were supposed to live in 
this state until disturbed by a "merge r". Indeed, merger s are very 
energetic events that can shock heat ([Bu rns et al. 1997) and mix 
the ICM dGomez et al.ll2002t iRitchie & Tho mas 2002): through 
these processes they were supposed to efficiently destroy cool- 
ing flows. After the mergers, clusters were supposed to relax and 
go back to the cooling flow state in a sort of cyclical evolu- 
tion. With the fall of the "cooling flow" brough t about by the 
XMM-Ne wton and Chandra observa tions (e.g. iPeterson et al.l 
1200 ll and lMolendi~& Pizzolato 2001), doubts were cast also on 
the interpretation of mergers as the dominant mechanism which 
could transform CC clusters into NCC. More generally speaking, 
the question arose whether the observed distribution of clusters 
was due to a primordial division into the two classes rather than 
to evolutionary differences d uring the history of the clusters. 
iMcCarthv et al.ld2004 12008) noticed the absence of systems that 
resemble observed NCC clusters in cosmological simulations, 
suggesting that mergers could not be the origin of the cluster dis- 
tribution. They proposed that early episodes of non gravitational 
pre-heating may explain the dichotomy; they envisage a scenario 
in which NCC clusters have been pre-heated to levels greater 
than ~ 300 keV cm 2 and do not have enough time to develop a 
cool core while CC clusters have been pre-heated to lower lev- 
els and need an additional source of present-day heating to off- 
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set cooling. With this model. [McCarthy et al.l d20 08) explained 
observed entropy and gas density Cha ndra profiles for CC and 
ROSAT gas density profiles for NCC. lQ'Hara et al.l (120061) sup- 
ported the primordial model showing that the scatter in scaling 
relations is larger for CC clusters than for NCC, suggesting that 
CC are not more relaxed than NCC system s . Mo reover the two 
body idealized simulations by iPoole et ail d2008l) showed that 
mergers cannot produce extended "warm" cores and cannot de- 
stroy metal abundance gradients. 

However, the evolutionary "merger" scenario ha s been continu- 
ously supported by observations. For instance, [Sanderson et al.l 
( 2006) compared temperature and cooling time profiles of a sam- 
ple of clusters observed by Chandra with indicators of merger 
activity and suggested that merger may be the pr imary factor in 
preven ting the formation of CC. More recently. ISanderson et alj 
(2009a|) have shown in a large sample of objects that the X- 
ray/BCG projected offset correlates with the gas density profile, 
which can be considered an indicator of the CC state. If the X- 
ray/BCG offset "measures" the dynamical state of the cluster, 
this result implies that the cool core strengths diminishes in more 
dynamically disturbed clusters. Another indicator of dynamical 
activity is the presence of extended radio halos on Mpc scales, 
whose ori gin is likely related to turbulent acceleration driven by 
mergers (Ferr ari et al.l 120081 for a recent review). None of the 
clusters which are found to host a radi o halo (on Mpc scale s) in 
a complete survey (IVenturi et al.l l2008: Bru netti et al.ll2.009l) can 
be classified as a CC object. These relations are indeed expected 
if mergers can efficiently destroy cool cores. 
An interesting model has been proposed by iMotl et all (|2004), 
who suggested a double role for mergers in the cluster formation 
process: while they can shock heat the ICM they also efficiently 
mix the ICM and participate to the formation of CC by provid- 
ing cool gas. By analyzing simulated temperature and surface 
brightness maps, they find that observational signatures of a cool 
core may disappear after mergers but cool gas remains in the 
systems at all time. Indeed mergers are complicated phenomena 
during which the ICM o f the in teracting objects is mixed. Also 
Ascasibar & Marke yitohl d2006l) . with different simulations with 
the purpose of studying the formation of cold fronts, showed that 
an interacting subcluster may donate its cool gas to the main 
cluster and that, during the merging processes, many hydrody- 
namic effects contribute to the mixing of the gas. 
Recently, we found an unexpected and interesting result from the 
analysis of the metalicity profiles in a large sample of clusters 
(Paper I): some clusters that cannot be classified as CC show a 
metal abundance excess at their center, witout a significant tem- 
perature decline or an X-ray brightness excess. We suggested 
that at least some NCC clusters have spent part of their lives as 
CC objects and therefore there cannot be a complete primordial 
separation of the two classes of objects. 

However, further analysis is needed on this result to use it to 
gain insight on the origin of the distribution of CC-NCC ob- 
jects. We would like to better characterize these regions and to 
know if and to what extent they are common in the population 
of NCC objects. The results in Paper I, as well as those de- 
scribed above comparing observations and simulation, are based 
on "uni-dimensional" properties (entropy and metal abundance 
profiles). Indeed, global properties and unidimensional profiles 
are easy to derive and it is possible to compare them quantita- 
tively but they inevitably result in a loss of information. There 
is now a large number of thermodynamic maps in the literature, 
which show that spherical symmetry is not generally fulfilled 
in clusters and revealed the presence of region s with "anoma- 
lous" characteristics outside of the cores (e.g lSun etai]l2002t 



iRossetti et al.ll2007l and many others ). 

In this paper, we present a systematic two-dimensional analysis 
on a sample of 35 clusters observed with XMM-Newton, to look 
for and characterize regions with a significant metal excess in 
NCC clusters, as those found in Paper I. The outline of the paper 
is as follows: in Sect. 2 we describe the sample and the data anal- 
ysis, in Sect. 3 we discuss the classification schemes and we de- 
fine "CC remnants" in Sect. 4. We interpret our results in Sect. 5 
and we summarize our findings in Sect. 6. Quoted confidence 
intervals are 68% for one interesting parameter unless otherwise 
stated. All results are given assuming a ACDM cosmology with 
Q m = 0.3, Q A = 0.7, and H Q = 70 km s _1 Mpc" 1 . 

2. Data analysis 

2.1. The sample 

The starting point of the sample of galaxy clusters d escribed in 
this p aper is the "B55" X-ray flux limited sample (Edg e et ail 
1990). In order to have a significant coverage of the cluster 
within EPIC field of view, we have eliminated the nearest 
clusters and considered only those with redshift z > 0.03. Then 
we analyzed all public XMM-Newton observations, including 
mosaics and multiple observations, as described in Sect. 12.21 
We discarded observations where, after soft proton cleaning, 
the effective exposure time (MOS 1 + MOS2 + pn) is lower than 
25 ks, except in the case of mosaics where observations with 
t exp < 25 ks have been used for images but not for spectral 
analysis. 

The list of the clusters which survived our selection criteria 
is given in Table Q] A644, A2244 have not been observed by 
XMM-Newton, while all the observations of A3391, A1736, 
A2142, A2063 and A1651 are badly contaminated by soft 
protons (t exp < 25 ks) and have been discarded. The clusters 
Cygnus A has been discarded in a second phase of the analysis 
because of the presence in the core of the radio galaxy QSO 
B1957 + 405, featuring large hotspots well detected in X-rays, 
which cannot be easily subtracted from the core spectrum. 



2.2. General data reduction 

We retrieve Observation Data Files (ODF) from the XMM- 
Newton archive and process them with SAS software version 
7.0. The event files produced by this standard analysis technique 
are then cleaned to remove soft proton flares with a double fil- 
tering process. As a first step, we produce the light curve in a 
hard energy band (10-12 keV) and we remove all time periods 
with count rates exceeding a fixed threshold 0.025 cts/s for the 
MOS detectors and 0.050 cts/s for the pn. This procedure al- 
lows to remove most flares, but softer flares may survive. In the 
second step, we apply a cr clipping technique to the histogram 
obtained from the light curve in 2 - 5 keV energy range. For 
each observation, we have c alculated the "in over out ratio", Rsb 
(iDe Luca & Molendlll2004l) . to identify observations badly con- 
taminated by a quiescent soft proton component. As outlined in 
Paper I, the ratio has been calculated in an external annulus at 
E > 9 keV to reduce the contribution of the cluster emission 
which fills all the FOV, especially for the nearest objects. We 
have discarded only one observation (and therefore the cluster 
A2065) where the R S b of the two MOS is larger than 2.0. After 
soft proton cleaning we filter the event files, according to pattern 
and flag criteria, and we remove by eye brightest point sources. 
Since we are mainly interested in characterizing the thermody- 
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namic properties of the ICM in the central and brightest regions 
of the cluste rs of our sample, advanced procedures to treat the 
background dLeccardi & Molendill2008bl) are not strictly neces- 
sary. This is also the reason why we could discard only badly 
contaminated observations with Rsb > 2.0. As background event 
files, we merged nine "blank sky" field observations, as com- 
monly done, and we extract images and spectra for the back- 
ground as for the source observations. We calculate a normaliza- 
tion factor Q for each cluster observation, to take into account 
possible temporal variations of the instrumental background. 
The normalization factor is the count rate ratio between source 
and background observations in an external ring (10' - 12') be- 
yond 9 keV. Background images and spectra are scaled by Q 
before subtraction from source images and spectra. 

2.3. Two-dimensional analysis 

For all the clusters in our sample we prepared two dimensional 
maps of the main thermodynamic quantities, starting from 
EPIC images. To do this we have used a modified version 
of the ad aptive binning + b road band fitting technique de- 
scribed in Rossetti et al. ( 2007t) . where we have substituted the 
Cappellari & Copini (l2003h adap tive binning algorithm with the 
weigh ted Voronoi tessellation bv lDiehl & S tatler 20061 dRossettil 
120061) . 

Especially in clusters undergoing major mergers, where there is 
no spherical symmetry, maps are a fundamental tool to select 
interesting regions for a proper spectral analysis, which is 
necessary to complement the thermodynamic al information 
with chemical information. Indeed a "blind" spectral analysis in 
concentric annuli would not allow to detect interesting features 
in our clusters. 



2.4. One-dimensional analysis 
2.4.1. Profiles 

As a first step, we have performed spectral extraction and anal- 
ysis in concentric annuli. As discussed in Sect. I2.3I the main 
drawback of this technique is that the assumption of spherical 
symmetry is not always fulfilled in our clusters and in some 
cases even the choice of the center can have a strong impact 
on the observed properties. Therefore we have used the thermo- 
dynamic maps (Sect. \2.3\ to identify those clusters where the 
deviations from the spherical symmetry are larger and we have 
decided not to perform radial a nalysis in four well known merg- 
ing o bjects: A3667 dBriel et al.ll2004l). A2256 dSun et al Il2 002), 
A754 (Henry et al.ll2004l) and A3266 dFinoguenov et al.ll2006l) 
For the remaining clusters of the sample, where we do not ob- 
serve large displacements between the surface brightness peak 
and the entropy minimum, we select as center of symmetry the 
surface brightness peak, even if the large scale isophotes have 
another center. This allows to better describe the ICM properties 
in the more central regions of the clusters. 
We extract spectra from annular regions around the selected cen- 
ter. For each instrument (MOS1, MOS2 and pn) and each region 
we extract source and background spectra and we generate an 
effective area (ARF) file. Then we associate to the spectrum a 
redistribution matrix file, appropriated for the instrument and, in 
the pn case, for the position of the selected region in the detec- 
tor. 

We perform spectral fitting, using the XSPEC vll.3 package, 
for each spectrum separately in the energy range 0.5-10 keV. 



We use an absorbed mekal model (WABS*MEKAL), where the 
Nh is fixed to the input galactic valu^B and the redshift is al- 
lowed to vary in a small range (width ^ 0.02) around the nom- 
inal value. Temperature, metal abundanc^ and normalization 
are the free parameters of t he fit (following the prescription in 
lLeccardi"& Molendi 2008a even negative values are allowed for 
the metal abundance). Best fit results obtained from the three 
instruments and from multiple observations of the same region 
are then averaged together with a weighted mean. This allows 
to produce projected temperature, metal abundance and surface 
brightness profiles for each cluster. 

Finally, we perform deproj ection of the pr ofiles, following 
the technique described in lEttori et all d2002l) . to derive three- 
dimensional density and temperature profiles. Combining them, 
we obtain entropy and cooling time profiles. 

2.4.2. IN and OUT regions 

In order to classify clusters according to the scheme described 
in Paper I, we extracted spectra in a IN region and in a OUT 
reference region, whose radii are defined as a fixed fraction of 
Rm (r < 0.05 Rm and 0.05^igo < r < 0.2 R m , respectively). 
R\so has been calculated as 



R 



180 



1780 



kT \ 
5keV/ 



1/2 



h(z) 1 kpc, 



(1) 



where h(z) = (Q m (l + z) 3 + Qa) 1 / 2 and kT is the temperatur e 
of the cluster dArnaud et alj|2005t iLeccardi & Molendil l2008b). 
This mean temperature has been first calculated starting from the 
profiles (eventually excluding the bins of the temperature drop 
in CC clusters). After extraction of the OUT spectra, we have 
re-calculated R^o using as mean temperature T out , i.e. the tem- 
perature obtained with a spectral fit in the region 0.05 R\w < 
r < 0.2/?i8(). We have then calculated the differences between 
the new estimate of 0.2R\%q and the old one. Once converted in 
arcseconds, we have compared this difference with the XMM- 
Newton point spread function: if the two estimates of Q.2R\%q 
differed by more than 15", we extracted again spectra with the 
new /?i8o and iteratively repeated the procedure until differences 
were smaller than the PSF In the case of the four clusters with- 
out radial profiles (A754, A2256, A3266 and A3667), the first 
estimate of the mean temperatures has been performed directly 
from the temperature maps. The list of the final T out and/?igo can 
be found in Table [2] 

As discussed in Sect. 12.31 the choice of the center has been per- 
formed starting from the thermodynamic maps. More specifi- 
cally, we have selected as a center the position of the minimum 
in the pseudo-entropy map. In most cases, the position of the 
entropy minimum coincides with the surface brightness peak, 
while in other clusters there is a significant displacement, but 
usually smaller then the radius of the IN regions, with the excep- 
tion of A3667. Moreover the choice of the position of the entropy 
minimum as a center instead of the surface brightness peak does 
not alter significantly the value of the pseudo entropy ratio, ex- 
cept for the case of A3667. This is a wel l studied merging cluster 
dVikhlinin et al.ll200U feriel et al.ll2004l) : the low-entropy gas is 
concentrated at the position of a prominent cold front, about 500 
kpc SE from the surface brightness peak (see Figure in Appendix 
B). 



1 The only exception is A478, where we found a difference of about 
a factor of two between the galactic value (1.48 1 21 cm 2 , Dickey & 
Lockman) and the best fit value (2.7 10 21 cm 2 ) 

2 For solar abundances we refer to Anders & Grevesse ( 1989) 
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Figure 1 : Comparison of temperature and emission measure ra- 
tios for the clusters in our sample. The solid line represent the 
threshold used to divide clusters between low entropy clusters 
and non low entropy clusters, corresponding to sin/sout = 0.45. 
The symbols and colors are as in Table Q] 

3. Classification schemes 

As in Paper I we plot the temperature ratios (T IN /T 0UT ) versus 
the emission measure ratios (EMm/EMout) f° r the clusters of 
our sample (Fig.Q} and we use pseudo-entropy ratios to divide 
clusters into two classes (Table |5J. We recall here the definition 
of the pseudo-entropy raticU 

^^^W" 3 (2) 
sout Tout \EMout ) 

The pseudo-entropy ratio is well correlated with the entropy 
ratio (see Appendix A) and therefore it is a useful and easy- 
to-calculate indicator of the variation of the physical three- 
dimensional entropy in the cores of galaxy clusters. 
As in Paper I, we compare our pseudo-entropy classification 
with two alternative classification schemes, based on the core 
properties (Col. 6 in Table |2]l and on the dynamical state (Col. 
7 in Table 0. More specifically, in Col. 6 we divide clusters in 
three classes: cool core (CC), intermediate systems (INT) and 
non cool core (NCC), where CC feature a prominent surface 
brightness peak and a temperature gradient, NCC possess nei- 
ther of these properties, while INT show only one of these ob- 
servational features. In Col. 7 we divide clusters in two classes: 
major mergers (MRG) and clusters showing no evidence of a 
major mergers (NOM). We consider evidence of a major merger 
the presence of cluster- wide diffuse radio emission, multi peaked 
velocity distribution of galaxies and significant irregularities on 
the surface brightness and temperature maps. The lack of these 
properties is not sufficient to state that an object is relaxed, this is 
why we prefer to consider these objects as "no observed merg- 
ers". We refer to Paper I for more details on the classification 

3 X-ray astronomers usually define "entropy" the quantity S = 
Txite 213 , where n e is the electronic densit y and T x the deprojected tem- 
perature of the ICM dPonman et al.ll999) . In this paper, we use the def- 
inition of "pseudo-entropy" s = T * (EM) ", where T is the projected 
temperature and EM is the emission measure, i. e. t he normalization of 
the MEKAL model in XSPEC per units of area (see lRossetti et al.l2007l 
for more details). 



schemes and for the necessary references. 
For the purposes of this paper, we are interested in separating 
clusters with a low entropy core (LEC) from the rest, this is why 
we have defined only one threshold to divide LEC from non- 
LEC objects. With our selected threshold (sjn I sout = 0.45), we 
can state that all LEC clusters are known to be CC and do not 
show any evidence of merger (except A85, probably undergoing 
a merger in an early stage, which has not affected the observa- 
tional properties of the core, see Paper I). Concerning non-LEC 
clusters, they are almost all non cool core systems, and many of 
them show significant indications of a major merger. As already 
discussed, the fact that some of them do not show significant 
merger features does not mean that they are relaxed. 

4. Cool core remnants 

The central regions of LEC clusters in our sample are charac- 
terized by low-entropy (by definition), often accompanied by 
a temperature decrease and a surface brightness excess. These 
regions are also characterized by a high metal abundance: the 
mean iron abundance in the cores of LEC clusters is signifi- 
cantly larger than the typical value of outer reg ions of galaxy 
clusters Z = 0.23 Zrr> dLeccardi & Molendi|l2008al) . even if with 
a a 20% scatter dDe Grandi & Molendil 120091) . It has been 
shown that the excess abundance is consistent with being due 
to metals ejected f rom the BCG ga l axy th at is invariably found 
in these systems dDe Grandi et alj 120041) . At present, we are 
not aware of any galaxy cluster with central regions charac- 
terized by low entropy gas without a metal abundance excess. 
Therefore, there must be a close relation between entropy and 
metal abundance, and generally speaking, between the thermo- 
dynamics and the chemical properties of the ICM. An exam- 
ple of this relation can be seen in Fig. [2] where, for each LEC 
cluster, we plot the metal abundance radial profile as a func- 
tion of the pseudo-entropy ratio. The pseudo entropy ratio is 
defined as in Equation |2j with the one difference that the tem- 
perature and emission measure of a given annulus replace Tin 
and EMof, i. e. s(r)/s 0U T = T(r)/T 0UT * (EM(r)l EM 0UT y l l\ 
In the lower panel of Fig.|2l we plot the mean error weighted 
abundance profile as a function of the pseudo-entropy ratio for 
LEC clusters and the one cr scatter around the average of the 
values. Low entropy ICM is characterized by a a significant 
metal abundance excess, although the scatter is quite large for 
*/ Sout < 0.4. The plots show that all LEC clusters show a sig- 
nificant metal abun dance excess with respect to the outer mean 
value Z — 0.23Zq dLeccardi & Molendi 2008a) in regions char- 
acterized by a pseudo entropy ratio smaller than 0.8 (correspond- 
ing to physical radii r < O.OlRiso)- 

In Paper I, we found that some clusters without a low-entropy 
core presented an unusually high metal abundance in their IN 
region (as already discussed such high abundances are typical of 
LEC). We concluded that the most likely explanation for the high 
central abundance of these anomalous non LEC systems was that 
at some time in the past they hosted a cool core (i.e. low-entropy 
gas and metal abundance excess) that was subsequently heated 
up. Indeed, a heating event that does not completely disrupt a 
cool core will most likely leave behind a region characterized 
by a high metalicity and by an entropy that, albeit not as low 
as the one found in the central regions of LEC systems, will be 
lower than that found in other regions of the cluster. We empha- 
size that metals are reliable markers of the ICM, in the sense 
that, once they have polluted a given region of a cluster, the time 
scale over which th ey will diffuse is comparable to the Hubble 
time (Sarazin 1988). Thus metals trace the ICM where they have 
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Table 2: Properties of the clusters in the sample: redshift (as derived from NED), "mean" temperature (i. e. Tout) and /?iso- Col. 5 
lists the values of sm/sout we have derived, Col. 6 and Col. 7 list alternative classification schemes (see text). The first 14 entries 
are LEC clusters, while the remaining are non-LEC. 



Name 


Redshift 


Temperature 
(keV) 


^180 

(Mpc) 


a 


Core CI. 


Dyn. CI. 


A2199 


0.0301 


4.101 


± 0.025 


1.58 


0.4216 i 


t 0.0033 


cc 


NOM 


A496 


0.0329 


4.990 


±0.038 


1.75 


0.3109 = 


t 0.0028 


cc 


NOM* 


2A0335+096 


0.035 


3.649 


± 0.009 


1.49 


0.2644 = 


t 0.0007 


cc 


NOM 


A2052 


0.035 


2.884 


±0.014 


1.33 


0.4274 = 


t 0.0028 


cc 


NOM 


A4059 


0.047 


3.968 


± 0.028 


1.55 


0.4394 = 


t 0.0044 


cc 


NOM 


Hydra A 


0.054 


3.413 


± 0.024 


1.43 


0.4138 = 


t 0.0040 


cc 


NOM 


A85 


0.055 


5.491 


± 0.048 


1.82 


0.3857 = 


t 0.0050 


cc 


MRG 


A1795 


0.062 


5.416 


± 0.026 


1.80 


0.3524 = 


t 0.0023 


cc 


NOM 


A3 112 


0.075 


4.319 


± 0.033 


1.59 


0.3855 = 


t 0.0038 


cc 


NOM 


A2597 


0.085 


3.491 


±0.018 


1.43 


0.3568 = 


t 0.0026 


cc 


NOM 


A478 


0.088 


7.057 


± 0.034 


2.03 


0.3706 = 


t 0.0024 


cc 


NOM* 


PKS0745-191 


0.103 


7.765 


i ±0.12 


2.11 


0.2957 = 


t 0.0056 


cc 


NOM* 


A2204 


0.152 


7.412 


± 0.092 


2.01 


0.2972 = 


t 0.0046 


cc 


NOM 


A2029 


0.077 


6.173 


± 0.056 


1.90 


0.4505 = 


t 0.0057 


cc 


NOM 


A4038 


0.030 


2.957 


±0.015 


1.35 


0.5093 = 


t 0.0037 


INT 


NOM 


A576 


0.039 


3.777 


± 0.044 


1.52 


0.6226 = 


t 0.0135 


INT 


MRG 


A3571 


0.039 


6.253 


± 0.045 


1.95 


0.6220 = 


t 0.0078 


INT 


NOM 


A119 


0.044 


6.032 


± 0.082 


1.92 


0.7621 = 


t 0.0279 


NCC 


MRG 


MKW3s 


0.045 


3.325 


± 0.018 


1.42 


0.4860 = 


t 0.0039 


CC 


NOM 


A 1644 


0.047 


4.156 


± 0.056 


1.59 


0.5318 = 


t 0.0132 


INT 


NOM 


A3558 


0.048 


5.221 


± 0.025 


1.78 


0.6155 = 


t 0.0057 


INT 


NOM 


A3562 


0.048 


4.344 


± 0.032 


1.62 


0.6089 = 


t 0.0083 


NCC 


MRG 


Triangulum Australis 


0.051 


9.224 


± 0.095 


2.36 


0.6845 = 


t 0.0164 


NCC 


NOM 


A3158 


0.060 


4.903 


± 0.042 


1.71 


0.7488 = 


t 0.0153 


NCC 


MRG 


A399 


0.072 


5.984 


±0.109 


1.88 


0.7086 = 


t 0.0319 


NCC 


MRG 


A401 


0.074 


7.268 


± 0.084 


2.07 


0.7028 = 


t 0.0212 


NCC 


MRG 


A2255 


0.081 


6.197 


±0.154 


1.91 


0.8853 = 


t 0.0714 


NCC 


MRG 


A 1650 


0.084 


5.427 : 


± 0.0539 


1.78 


0.5323 = 


t 0.0089 


INT 


NOM 


A 1689 


0.183 


8.612 


±0.102 


2.13 


0.4724 = 


t 0.0089 


INT 


NOM 


A2319 


0.056 


8.962 


±0.091 


2.32 


0.5643 = 


t 0.01 13 


NCC 


MRG 


A3532 


0.055 


4.839 


±0.104 


1.70 


0.7285 = 


t 0.0336 


INT 


NOM* 


A3667 


0.056 


5.461 


± 0.029 


1.81 


0.6378 = 


t 0.0078 


NCC 


MRG* 


A754 


0.054 


8.544 


±0.102 


2.27 


0.5289 = 


t 0.0120 


NCC 


MRG 


A3266 


0.055 


7.351 


±0.081 


2.10 


0.6772 = 


t 0.0170 


NCC 


MRG* 


A2256 


0.057 


5.709 : 


± 0.0763 


1.85 


0.6300 = 


t 0.0197 


NCC 


MRG* 



Not es * Properties and the references for the clusters not included in the sample in Paper I. A 496: CC and no evide nce for a major merger 
dTamura et alJl200H and references therein). A478: C C and no evidenc e for a major merger (Sand erson et alj [2005 and references therein). 
PKS0745-191: CC and no evidence for a major merger ( Chen et al. 2003 and references there in). A3532: INT (moderate temperature drop) and 
no evidence for a major merger (the indi cations of optical subst ructures, fBardelli et alj2000j are insuffic ient to claim a major merger). A3667: 
NCC and MRG ( evidence of radio reli cs Rottger ing et alj 19971 . optical substructures lOwers et alj, 2009. and disturbed ther modynamica l maps 
iBriel etalj l2004; Vik hlinin etaUl2001ri . A3266: NCC and M RG (even if extended radio emission h as not been detected. lBuotdl200ll . there 
are evidences of optical substru ctures. lOuintana et all H 996. and disturbed thermodyna mical maps iFinoguenov e"taUl2006h. A2256: NCC 
and MRG (radio halo and relic. I Kiml 1 19991 . optical substructures, iBerrington et al]|2002l . and disturbed temperature maps. [Sun et al.|[2002l ; 
iBourdin & MazzottiJl2008h . 



been injected. The presence of these regions h as also been pre- 
dicted by the simulations of Mot l et all d2004l) . who suggested 
that even if the observational cool core signature may disappear, 
cool gas (and we may add "low-entropy and metal rich gas") re- 
mains in the merging systems at all times. 
Since in LEC clusters a metal abundance excess is invariably 
associated to ICM with a low pseudo-entropy, we have system- 
atically selected the regions with the lowest pseudo-entropy ra- 
tio in non-LEC systems, with the aim of finding regions pos- 
sibly characterized by a high metalicity. In practice, for the 
21 non-LEC clusters of our sample, we have prepared two- 
dimensional pseudo-entropy ratio maps, by dividing the pseudo- 
entropy maps for the pseudo-entropy in the OUT region derived 
with spectral analysis. In these maps we have identified regions 
characterized by a s/squt < 0.8 (as shown in Fig. [2] all LEC 



clusters show a metal abundance excess in regions with entropy 
ratio < 0.8), and we have extracted and analyzed spectra in these 
regions (as described in Sect. 12.4. U to investigate their metalic- 
ity. Pseudo-entropy ratio maps and figures of the selected regions 
are provided in Appendix B. 

In all the objects (21) of our subsample of non-LEC clusters, we 
find regions characterized by an entropy ratio smaller than 0.8. 
In Fig. [3j we plot the metal abundance in these regions for all 
the clusters (results are reported in TableQ. 12 objects (namely: 
A1689, MKW3s, A4038, A1650, A1644, A3558, A576, A754, 
A3562, A3571, A3667 and A2256) show a significant excess 
(> 3cr c.l.) with respect to the reference value Z = 0.23Zq 
dLeccardi & Molendill20 08a). In the remaining 9 clusters, the 
metal abundance is statistically consistent with the reference 
value. Amongst these objects, we will consider as clusters with- 
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Figure 2: Upper panel: Metal abundance profiles of LEC clus- 
ters as a function of the pseudo-entropy ratio, in the regions of 
spectral extraction. The symbols and colors are as in Table Q] 
Lower panel: Mean error weighted abundance profile for LEC 
clusters. Error bars show the (small) error on the average, while 
dotted lines show the lcr scatter around the mean. In both pan- 
els, the horizontal dashed line indicates the mean value of outer 
region s of galaxy clusters Z = 0.23Z© dLeccardi & MolendU 
2008a) and the vertical dashed line indicates the selected thresh- 
old of s/squt under which all LEC cluster show a significant 
metal abundance excess. 



out metal abundance excess those where the metalicity is mea- 
sured with a good statistical accuracy, more specifically the 6 
clusters where Z > QAZq can be excluded at more than 2cr con- 
fidence level (A2319, A3266, A399, A401, A119 and A3158). 
We will no more consider in this paper the three objects (A3532, 
Triangulum Australis and A2255) where the error bars on the 
metal abundance are too large to discriminate between the two 
classes (the observations of these clusters are the shortest of our 
sample, see Table 1). 

The most likely interpretation is that the low entropy high Z re- 
gions found in many of our non-LEC systems are what remains 
of a cool core after a heating event. We dub these structures "cool 



Figure 3: Metal abundance in the regions selected for spec- 
tral analysis for the 21 non-LEC objects in the sample (upper 
panel) and as a function of the pseudo-entropy ratio in the re- 
gions (lower panel). The horizontal line represents the reference 
value for the metalicity of of outer reg ions of galaxy clusters 
Z = O.23Z (iLeccardi & MolendU l2008al) . Symbols and colors 
are as in Table Q] 



core remnants'Q a more detailed discussion of how they formed 
and of alternative scenarios is provided in Section 
In the lower panel of Fig. [3] we plot the metal abundance as 
a function of the pseudo entropy ratio, s/squt where s is the 
pseudo-entropy in the region used for spectral analysis. It is 
worth noting that the clusters with no indication of a signifi- 
cant metal abundance excess are those where the selected re- 
gions show the highest s/sout, with the notable exception of 
A2319 (red filled square) 

To better characterize the systems in our sample, we have com- 
pared the entropy maps with the galaxy distribution and more 
specifically the low-entropy regions with the position of the 
BCG (Col. 3 in Table [3] and Figures in Appendix B). The 
name and positions of the BCGs are derived from the literature 
(ICoziol et al.ll2009t iLin & Mohrll2004t iGudehus & Hegvill 1991b 
iPostman & Lauerlll995l) . In 9/12 clusters with "CC remnants" 



4 The term "remnant of cool cores" has been recently used by 
iFusco-Femiano e t al. ( 2009) to describe structures similar to the ones 
we detect albeit without information on the metal abundance. 
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(namely A1650, MKW3s, A4038, A1644, A3558, A3562, 
A3571, A 1689 and A576) the BCGs are located in the low- 
entropy high metal abundance regions. In A3667 and A2256, 
the positions of the BCGs do not coincide with the "CC rem- 
nant", which is moreover not obviously associated to any galaxy 
concentration (a more detailed s tudy of the galaxy distribution 
in A3667 by lOwers et al.ll2009l found essentially the same re- 
sults). In A754 the BCG is located ~ 700 kpc west of the se- 
lected regions but in the "CC remnant" we find the elliptical 
galaxy 2MASSX J09091923-0941591 (see Figure in Appendix 
B), which is the bri ghtest member o f the secondary peak in the 
galaxy distribution (Fab ricant et al.llT986b . The position of this 
secondary galaxy concentration, possibly associated to the CC 
remnant, coincides also with a clump in the dark matter distri- 
bution , as derived with gravitational lensing (lOkabe & Um etsu 
2008). For the non-LEC objects where we do not observe a sig- 
nificant metalicity excess, we find that in 3/6 systems (A119, 
A401 and A3 158) the low entropy regions are not obviously as- 
sociated to the BCGs. In the remaining 3 clusters (A399, A3266, 
and A2319) the BCGs are located in the low-entropy regions 
where we do not find evidence of a metal excess. 
The main reason for which we compared the position of CC 
remnants with the galaxy distribution is that galaxies, and more 
specifically BCGs, are respon sible of the metal abun dance ex- 
cess observed in LEC clusters (De Grandi et al . 2004). However 
the galaxy distribution may also provide information on the 
shape of the potential well of the clusters. Indeed, in relaxed 
clusters BCGs are located at the bottom of the potential well and 
their positions coincide with the centroid of the X-ray isophotes 
and with the brightness peak of their host clusters. We have com- 
pared the position of CC remnants with the position of the X- 
ray peak, the centroid of X-ray isophotes and the position of the 
BCGs (Figures in Appendix B). In most cases (namely A1650, 
MKW3s, A4038, A1644, A1689, A576, A3562, and A3571), 
these three points lie within the low entropy regions even if they 
do not necessarily coincide with one another. However, in the re- 
maining clusters one or more of these indicators is significantly 
offset from the CC remnant. A754 and A2256 have a very dis- 
turbed morphology and the X-ray centroid does not coincide 
with the CC remnant or with the BCG. In A3667 the centroid 
roughly coincides with the brightness "peak" and the BCG, but 
the CC remnant is offset by almost 500 kpc. It is apparent that in 
these last three cases, the low entropy metal rich regions are not 
in equilibrium within the cluster's potential well and therefore 
we are observing a rapidly evolving situation. 

5. Discussion 

In Section|U we have shown that 12 out of 21 non-LEC objects 
of our sample, host regions characterized by low entropy and 
high metal abundance, with respect to mean values in the outer 
ICM. Such conditions are usually verified in the cores of relaxed 
LEC clusters, albeit with stronger gradients (i.e. the entropy of 
the ICM reaches smaller values). We have dubbed these features 
"cool core remnants", since we have interpreted them as what 
remains of a cool core after a "heating event". 
One may argue that these structures have evolved similarly 
but independentely of LEC systems and that there is no need 
to trace them back to known structures, such as LEC. In the 
"primordial" scenario, where non-LEC clusters started off with 
different early condition and are now slowly cooling to become 
LEC, low entropy and metal rich region could be interpreted as 
"progenitors" of cool cores, rather than remnants. However this 
scenario does not explain why in some cases these structures 



Table 3: Properties of low entropy regions. In col. 2 we list the 
metal abundance in solar units, while in col. 3 we write "Y" if 
the BCG is found in the low-entropy region and "N" if it is else- 
where. Clusters in italic are "AGN CC-remnant" while clusters 
in bold type are "merger CC-remnant" (Sect. [5] for these defini- 
tions). 



I 1 1 1 c t&r 


\A t^tai A r\i lnnQnpp 1 r tW or) 
IVlCLcll rtUUIlilallLC ^sUlcll ) 


D \_, VJ 


A 

r\l OJU 


O ^Q7 -i- O CO 1 
U.J7 / ± U.UZ1 


v 

I 


IVIA. WjS 


n A^Q -u- o ono 


V 


/i'+L/Jo 


U.^fUZ ± U.UUo 


v 


A 1AAA a 




v 

1 


A1007 


KJ.Zyo + u.ui / 


v 

1 


A 


n 4.nnin 


v 

I 


a zn(\ 

t\.j 1 o 


U.jZU ± U.UjJ 


1 


A/54 


r\ 1 1 a _i_ n no i 
U.314 ± U.UZ1 


Mb 


A3562 


0.471 +0.021 


Y 


A3571 


0.349 + 0.013 


Y 


A3667 


0.327 +0.012 


N 


A2256 


0.351 +0.038 


N 


A2319 


0.242 + 0.012 


Y 


A3266 


0.232 + 0.039 


Y 


A399 


0.221 +0.053 


Y 


A401 


0.262 + 0.044 


N 


A3158 


0.275 + 0.077 


N 


A119 


0.242 + 0.078 


N 


A2255 c 


0.547 + 0.224 




A3532 c 


0.473 + 0.095 




Triangulum Australis 


0.369 + 0.151 





Notes: a Even if the central entropy of this cluster is within reach 
of AGN outbursts, it also show some indications of on-going in- 
teraction and cannot be unambiguously associated to the "AGN 
CC-remnant class" (see text). b The BCG 2MASSX J09083238- 
0937470 is not associated to the CC-remnant, which however 
coincides with the position of the elliptical galaxy 2MASSX 
J0909 1923-094 1591. c Due to the large indetermination on the 
metal abundance estimate, we could not unambiguously detect or 
exclude the presence of a metal abundance excess in this object. 



are not in equilibrium within the potential well of their host 
clusters and are not associated to the BCG or to a giant elliptical 
galaxy. On the contrary, these issues are naturally addressed in 
the "evolutionary" scenario if we consider the effects of mergers 
both on the ICM and on the galaxy distribution (Sec. I5.21 >. 
Moreover we should pay attention to the fact that most of these 
structures are found in clusters undergoing major mergers, i. e. 
in rapidly evolving objects. It is therefore necessary to trace 
these transient structures back to an equilibrium state from 
which they recently evolved. 

Another possible objection is that the low-entropy structures 
may not be embedded in the ICM, but rather in groups along 
the line of sight of the clusters. However, in the cases where the 
cool core remnants are associated to the BCGs or to a galaxy 
concentration (see Sect. @), the redshifts of the galaxies are 
consistent with the mean readshift of the clusters. In the cases 
where there is no galaxy concentration obviously associated 
to the low entropy regions, these features would have to be 
interpreted as "gas-only" groups superposed on the line of sight. 
In the interpretation of the low-entropy regions as CC remnants, 
the "heating event" is responsible for smoothing the entropy 
gradient of the cores of LEC clusters. We have identified two 
possible mechanisms: interaction with the central AGN and 
mergers. 
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5.1. CC remnants and central AGNs 



Interactions between the central AGN and the ICM have been 
observed with Chandra and XMM-Newton in many LEC clus- 
ters and these interactions are now considered as the prin- 
cipal mechanism preven ting the formation of cooling flows 
dPeterson & Fabianl2006l and references therein). In some cases, 
powerful AGN outbursts may increase significantly the entropy 
content of some galaxy clusters, transforming a low entropy 
core into a non-LEC object. This is possible for systems where 
the core entropy has been raised to a rel atively small value 
so * (30 - 50)keVcm 2 dVoit & Donahue! 120050 . In our sub- 
sample ofnon ; LEC_cjusters, according to Chandra measure- 
ments dCavagnolo et all I2009Q . only A1650, A4038, MKW3s 
and A 1644 have core entropy smaller than 50 keV cm 2 and there- 
fore are within the reach of powerful AGNs. AGN heating is 
the most plausibl e explanation for incr e asing the entropy in the 
cores of A1650 dDonahue et all 120051: IVoit & Donahue! 120051) 
and for MKW3s (radio lob es have been observed in this clus- 
ter, iGiacintucci et alj[2006h . We will consider these clusters as 
"AGN CC-remnants" 

. The remaining two non-LEC clusters with a relatively low cen- 
tral core entropy, A4038 and A1644, are controversial objects, 
showing some indications of ongoing interactions, contrary to 
A1650 and MKW3s which look relaxed at all wavelengths (see 
Table 2 in Paper I). A4038 does not have a large temperature 
drop at the center but it has a low central co oling time. This is 
why some authors classify it as a cool core dPeres et all 11998b 
and o thers as a non cool c ore dSanderson et al.ll2006l) . In the op- 
tical, iBurgett et alj d2004l) found some indications of substruc- 
ture on large scales. A1644 is probably an advanced off-axis 
merger with two subclumps clearly visible in the X-ray image 
and some indicatio ns of on-going sloshi ng of the cool core of 
the main subclump (Reiprich et al. 2004). The interpretation of 
A4038 and A1644 is not straightforward, since merging events 
could have contributed to the heating of the cores (see Sect. 15. 21 ). 
In AGN CC remnants, the low entropy regions coincide with the 
centers of the clusters and with the positions of the BCGs as 
expected. In Fig. |4]we compare the error weighted mean metal 
abundance profile of LEC clusters with the same profile for non- 
LEC clusters with sq < 50keVcm 2 . In the innermost bins (low- 
est entropy ratios), the mean profile of this class of clusters is 
significantly larger than in LEC clusters, although the scatter is 
large. This is consistent with a scenario where a localized heat- 
ing in the inner regions of clusters has increased the entropy ratio 
without substantially modifying the metal abundance, resulting 
in ICM with enhanced metalicity for its entropy. 



5.2. CC remnants and mergers 

As shown in Sect. 15.11 8 out of 12 clusters featuring a metal 
abundance excess have core entropies too large to be pro- 
duced even by the most powerful AGN outbursts. Therefore, we 
need another mechanism which could have produced the CC- 
remnants that we observe today. The fact that most of these 
clusters (namely A754, A3667, A2256, A3562 and A576) show 
strong indications of a major merger suggests that merging could 
provide the necessary heating to produce these features. The 
indications are not as strong for the remaining three clusters 
(A1689, A3558 and A3571; in Table|2]they are classified as "in- 
termediate" and "no observed merger"). However, as shown in 
the notes reported below, they all possess some peculiar features, 
which suggest that they are undergoing some kind of interaction. 
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Figure 4: Error weighted metal abundance profile as a function 
of the pseudo entropy ratios for LEC clusters (blue squares) and 
for the four clusters with low central entropy (red circles). Dotted 
lines show the one cr scatter of the data. 



It is therefore natural to define this class of 8 objects as "merger 
CC remnants". 

Notes on A1689, A3558 andA3571 

A1689 This cluster has been considered for a long time by X- 
ray astronomers as an example of a relaxed cluster, due to 
its spherical shape (Figure in A ppendix B) and v ery peaked 
surface brightness profile (e. g. iPeres et al.ll998). However, 
analysis of the distribution of the velocities of the galax- 
ies showed tw o distinct peaks, su ggesting a line of sight 
superposition (Girardi & Mezzetti 200 l j). More recently, an 
analys is of XMM-Newton data by lAndersson & Ma deiski 
(2004) showed asymmetric temperature and redshift distri- 
bution which, combined with optical data, could suggest an 
on- going merger. 

A3558 Rosset tTet al.l (120071) classify this cluster as an interme- 
diate object, since it has some characteristics of CC but it 
also shows some indications of an on-going merger, like 
asymmetry in the temperature and entropy di stribution and 
subst ructures in the 3-D galaxy distributions dBardelli et al.l 
1998). Due to the very peculiar environment in which it is 
located (the core of the Shapley Supercluster), these features 
may be explained as the results of multiple minor mergers or 
of a past off-axis major merger with A3562. 

A3571 This cluster has a regular morphology in X-rays (Figure 
in Appendix B), but we do not observe a temperature de- 
crease in the center. T his is why it is often classified as a 
"non cool core" object dO'Hara et al. 2006; Sand erson et al.l 
2006) . On the bas is of its radio and optical properties, 
Ventu ri et al. I d2002l) suggest that this cluster is a late stage 
merger. 

We recall here the main results concerning the position of 
merger CC remnants with respect to the centroid of X-ray emis- 
sion and to the galaxy distribution (Sect. [4] and Table [3]). In 
5/8 clusters (namely A3558, A3562, A3571, A1689, A576) the 
BCGs and the X-ray centroid are located in the low-entropy high 
metal abundance regions. In A754 the CC remnant is associated 
to a secondary peak in the galaxy distribution (with a giant el- 
liptical galaxy), but it is also offset from the position of the large 
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scale X-ray centroid. In A3667 and A2256 the CC remnants are 
not obviously associated to any galaxy concentration and they 
are not found at the centers of their host clusters (see figures in 
Appendix B). These different results are not in contrast with our 
interpretation since the effect of the merging processes both on 
the collisional ICM and on the non-collisional dark matter and 
galaxies depends on many parameters such as the merger state, 
the mass-ratios and the impact parameter. Therefore, while in 
A3667 and A2256 the low entropy and high metal abundance gas 
may have completely decoupled from its galaxy concentration 
and possibly dark matter halo (such as in the "bullett cluster", 
IClowe et al.[|2006l) , in A3558, A3562, A3571, A1689, A576 and 
A754 this gas has not been completely displaced from its grav- 
itational potential well. Another possibility is that, in the latter 
systems, we are currently observing a late stage of a merger, 
where the low entropy high metal abundance gas is slowly mov- 
ing back to the center of the cluster. 

Two main mechanisms can contribute to the formation of a 
CC remnant during a merger event: shock heating and mix- 
ing. Mergers are expected to drive moderately supersonic shock 
waves in the ICM, which, being irreversible changes, increase 
the entropy of the system. If one of the merging subcluster had a 
CC before the merging event, shock heating may significantly in- 
crease the entropy of the core but in some cases this heating may 
not be sufficient to bring the entropy to the large values observed 
in some non-LEC clusters. With a very simplified calculation 
based on Rainkine-Hugoniot jump conditions, we estimate that 
shocks with Mach number 1-3 can increase the mean entropy in 
the cores only of a factor up to 1.8. Mergers also drive motions 
in the ICM in which low entropy metal-rich gas may be put in 
contact with high-entropy metal poor ICM. The mean entropy of 
the resulting mixed ICM will be lower than ambient entropy (but 
higher than typical core entropy) and its metal content could be 
larger than in the ambient gas. However, since there is a signifi- 
cant spread in the metal abundance values of cool cores, mixing 
may also efficiently disrupt small abundance gradients. We re- 
call here that metals are frozen in the ICM and therefore mixing 
with metal poor ICM is the only possible way of reducing the 
metal abundance. 

In Fig. |5]we compare the mean metal abundance profile for clus- 
ters with merger CC remnants with that of LEC clusters, simi- 
larly to Fig. [4]for AGN CC remnants. In this case the mean pro- 
files are consistent but the scatter in clusters with CC remnants 
is larger. Indeed, we did not expect the same results as in Fig. 
@] (i. e. gas with enhanced metalicity for its entropy) because, 
contrary to AGN heating, merger heating does not affect only 
the central regions of clusters and therefore the entropy ratios 
are modified in a non-trivial way. Moreover we cannot neglect 
the effect of gas mixing, both on the entropy and on the metal 
abundance. 

5.3. Non LEC clusters without metal abundance excess 

Of the 21 objects of our sample which are classified as non LEC, 
6 features regions with pseudo-entropy ratios smaller than 0.8 
but without a significant metal abundance excess (Z < 0.4Z© at 
more than 2cr confidence level). We will discuss in this section 
some of the properties of these systems. 

As shown in Table 12 we have strong indications of ongoing 
mergers for 5 out of 6 objects, namely A2319, A399, A401, 
A119 and A3266. The remaining cluster, A3158, shows some 
indications of ongoing interactions based mainly on the galaxy 
distribution dJohnston-Hollitt et alj|2008l) . These clusters do not 
have a CC remnant, and this could mean that they have never de- 
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Figure 5: Error weighted metal abundance profile as a function 
of the pseudo entropy ratios for LEC clusters (blue squares) and 
for the clusters with merger CC remnants (red circles). Dotted 
lines show the one cr scatter of the data. 

velope d a cool core, as in scenario proposed bv lMcCarfhv et al] 
(2008). However we should note that while in this "primordial" 
scenario, NCC may indifferently be relaxed or non relaxed ob- 
jects, in most objects of this class (if not all) we have found in- 
dications of on-going interactions. 

The fact that most of these clusters show indications of on-going 
interactions, could suggest that they may have had a low-entropy 
core that has been completely erased during the merger event. As 
shown in many simulation works, the effects of a merger on the 
structure of the ICM depends on many parameters, such as the 
mass ratio, the impact parameter and the structure of the interact- 
ing subclusters. Therefore it is not surprizing that in some cases 
the merger can efficiently destroy the low entropy core, while in 
others it leaves a CC remnant. 

We recall here than in 3 clusters (A399, A3266 and A2319) the 
BCGs are located in the low-entropy region where we do not 
find evidence of a metal excess (see Table[3]). The interpretation 
of this class of clusters is an interesting task: we speculate that 
either the LEC state from which they evolved was characterized 
by a metal abundance profile showing only a moderate excess 
or, alternatively, that the mixing may have been more effective in 
these clusters than in those where we observe CC remnants. We 
recall that metalicity profiles of LEC clusters show a large scat- 
ter: as can be seen also in Fig. [2] while the central metal abun- 
dance of some clusters can reach almost solar values, in others it 
reaches only 0.4 Zq. Assuming a constant degree of mixing in all 
cluster, it is clear that it would be easier to detect a metal excess 
in regions where gas with Z ~ Zq has mixed with the ambient 
gas, than in regions where mixing has occurred between ICM 
with Z =i 0.4 Zq and ICM with Z 0.2 Zq- Therefore systems 
like A399, A3266 and A2319 could be remnants of "metal poor 
CC", where "poor" means that the metal abundance is lower 
with respect to other CC objects. Alternatively, these could be 
systems where metals have been mixed more efficiently than in 
others where an abundance excess is detected. 
It is interesting to note that the fraction of objects where the BCG 
is not associated to the low-entropy regions is higher in clusters 
showing no Z excess (3/6) than in merger CC remnants (2/8), 
even if we should note that this difference is not statistically 
significant. As shown in many simulations and in the case of 
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Figure 6: Metal abundance profiles of non-LEC clusters as a 
function of the pseudo-entropy ratio, in the regions of spec- 
tral extraction. The symbols and colors are as in Table Q] The 
dashed horizontal line indicates th e mean value of outer region s 
of galaxy clusters Z = O.23Z , dLeccardi & Molendil l2008a). 
while the vertical dashed lines marks the threshold pseudo en- 
tropy ratio used to select regions for spectral analysis. 



the "bullett cluster' ' (IClowe et alJ l2006>. in violent major merg- 
ers the collisionless galaxy population may completely decouple 
from the collisional ICM. It is natural to assume that these vio- 
lent mergers are also more effective in mixing the gas and there- 
fore in completely erasing metal abundance gradients. In this 
scenario, the three clusters showing no metal abundance excess 
and no BCG associated, namely Al 19, A3158 and A401, should 
be those that have undergone the most destructive interactions. 
However this speculation needs further verification with a larger 
sample of clusters. 

5.4. Metal abundance and entropy distribution 

One may speculate whether our selection based on the pseudo- 
entropy ratio < 0.8 is effective in identifying all the regions with 
a large metalicity. Indeed, a possible counter-example to our 
choice of using an entropy ratio threshold for the selection of 
regions for spectral analysis would be the observation of a large 
metal abundance associated to high entropy gas and generally 
speaking the lack of anti-correlation between metal abundance 
and pseudo-entropy ratio. 

In Fig. [2] we have shown that the metal abundance excess in 
LEC clusters is typically found in regions with s/sout < 0.8. 
In a few cases, an excess is detected also in regions with 
s/sout ~ 1 but not at larger values. Moreover, the metal 
abundance is always anti-correlated with the entropy. 
We have performed the same test on the subsample of non-LEC 
clusters, using the results of spectral analysis in radial annuli 
centered on the surface brightness peak (Fig. |6). We do not 
include in this figures the four highly disturbed clusters for 
which we have not performed spectral analysis in radial annuli. 
It is interesting to note that we do not find a single region where 
a significant metal excess is associated to a pseudo entropy 
ratio larger than 1 . If a metal abundance excess is observed it is 
invariably associated to ICM with a low pseudo-entropy ratio. 



Figure 7: Cooling time in the inner 0.05 R\ go for the clusters of 
the sample for which we have performed radial analysis and de- 
projection. Horizontal dashed line is the Hubble time. 

5.5. Cooling times 

The results presented in this paper support the "evolutionary" 
scenario of the CC-NCC dichotomy: most non-LEC clusters 
have likely "evolved" from a LEC state. It is therefore natural 
to ask ourselves if this is a "one way trip", i. e. if once a clus- 
ter has undergone a merger it remains a non-LEC cluster for the 
rest of its life or it has the possibility of reforming a low-entropy 
core, as in the old cyclic vision of cooling flows-mergers. 
In Fig. Q we plot the cooling time in a central bin of radius 
0.05,Ri8o for the clusters of the initial sample for which we 
have performed radial analysis and de-projection (i. e. excluding 
A3266, A2256, A3667 and A754). The cooling time is gener- 
ally larger than the Hubble time for non-LEC clusters, including 
those of most clusters with CC-remnant (A576, A3571, A3562, 
A3558, we do not have this 3-D information for A3667, A754 
and A2256). Therefore these clusters would appear to be unable 
to develop a new low-entropy core. 

The cooling time shown in Fig. |7]is a mean value calculated in 
large region. If the mixing of the ICM is not completely effec- 
tive, clumps of cool gas may reside in these regions and in these 
clumps the cooling will be much more effective. In this case, 
we should consider the measured cooling time only as an upper 
limit and it is possible that the cooler gas could contribute to 
the formation of a low entropy core on time scales shorter than 
the Hubble time. Any future observation of these clumps of cool 
gas, will be useful to test the possibility of a "return journey" to 
the LEC state. 

6. Summary and Conclusions 

In this paper, we have presented a systematic two-dimensional 
analysis of the entropy and metal abundance of a sample of 
nearby clusters observed with XMM-Newton . 

- We have analyzed the observations of a sample of 35 clus- 
ters. Following the prescription in Paper I, we have classified 
14 of them as low-entropy cores. In the remaining 21 non- 
LEC objects we have performed a systematical analysis of 
the entropy maps and selected regions with pseudo-entropy 
ratio lower than 0.8, for a proper spectral analysis. All the 
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clusters in our subsample host regions with these character- 
istics. 

- In most non-LEC objects (12/21) the low entropy regions are 
characterized by a significant metal abundance excess with 
respect to outer regions of galaxy clusters. For 6 of the re- 
maining clusters, we can significantly exclude an abundance 
larger than QAZq, while for 3 objects the error bars are too 
large to discriminate. We have interpreted the low-entropy 
high metal abundance regions that we found in 12 clusters 
as what remains of a cool core after a heating event, and we 
have dubbed them "CC remnants". 

- For two clusters, A1650 and MKW3s, the most likely 
heating mechan is m is AGN feedback. As shown by 
IVoit & Donahuel d2005l) . clusters with core entropy < 
50keV cm 2 are within the reach of powerful AGN outbursts. 

- In 8/12 clusters with CC remnant, the core entropy is too 
large to be produced by giant AGN outbursts. In all these 
clusters we have found indication of on-going interactions 
and 5 of them are probably undergoing a major merger. 
Therefore the most likely heating mechanism for this class 
of objects is merging, which can shock heat and mix the 
ICM. In most of these clusters, the CC remnant is associ- 
ated to the BCG or to a giant elliptical galaxy which could 
be responsible for the production of the metal excess dur- 
ing the LEC phase. In A3667 and A2256 the merging may 
have completely decoupled the low-entropy metal rich gas 
from its associated galaxy concentration. In A754, A2256 
and A3667 the position of the CC remnant does not coincide 
with the centroid of the X-ray emission, indicating that the 
low entropy metal rich ICM is not in equilibrium within the 
potential well of the host clusters. 

- In 6 objects out of 21, we have not found a significant metal 
excess in the low entropy regions. Since most of these clus- 
ters show indications of major mergers, we have interpreted 
them as clusters where the metal abundance gradient has 
been erased by the mixing following mergers. 

The results presented in this paper strengthens those shown in 
Paper I. Indeed, in Paper I we found a metal abundance ex- 
cess in a small fraction of non-LEC clusters. Thanks to the two- 
dimensional analysis presented in this paper, we have found a 
metal abundance excess in most objects. We conclude that most 
non-LEC objects have spent part of their life as low-entropy 
cores. After a heating event the LEC signature disappears but 
the "CC-like" ICM (i. e. characterized by large metal abundance 
and low entropy, albeit not as low as in LEC) remains in the sys- 
tems and can be identified in two-dimensional entropy maps. In 
the framework of the alternative "primordial" scenario the low 
entropy metal rich regions could be interpreted as "progenitors" 
of cool cores, that are now evolving (slowly cooling and increas- 
ing their metalicity) to become LEC. However this scenario does 
not explain why most of these regions lie in dynamically ac- 
tive objects and at least three of them are found in a configu- 
ration of non-equilibrium, offset from the center of the clusters. 
Conversely, this issue is naturally addressed if they are the re- 
sults of a merger event in the evolutionary scenario. 
The results summarized above strongly support scenarios where 
cluster core properties are not fixed "ab initio" but evolve across 
cosmic time. 
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Figure A. 1 : Entropy ratio versus Pseudo-entropy ratio 



As expected, the entropy ratios are smaller than pseudo-entropy 
ratios, since projection effects smooth out gradients. 
In order to quantify the scatter around the best fit relation, we 
have applied to the residuals a maximum likelihood algorithm 
that postulates a par ent distribution describ ed by a mean and an 
intrinsic dispersion (Ma ccacaro et aljfl988h . To do this we have 
calculated an "effective error" on the y coordinate, which takes 
into account the errors on both coordinates, defined as: cr 2 e ^ - 

a 1 + b 2 cr\, where b is the best-fit value for the slope of the linear 
relation between the two data sets (in this case b = 0.792, Eq. 
lA.U . The mean value of the residual distribution is consistent 
with zero (0.021), with an intrinsic dispersion cr = 0.043. 



Appendix B: Figures of individual clusters 

Large figures, available at: 

http://www.iasf-milano.inaf.it/~rossetti/public/CCR/rossetti.pdf 
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Appendix A: Relation between Pseudo-Entropy and 
Entropy ratios 

We have used the entropy profiles for the 31 clusters in our 
sample for which we have performed radial analysis, to test 
whether the pseudo-entropy ratios (defined in Paper I) are good 
indicators of the behavior of the true three dimensional entropy 
(S = Tx/n^ 3 ). We have therefore calculated the entropy ratio 
S in /S out, starting from the three-dimensional entropy profiles 
obtained with the de-projected temperature and density profiles. 
We show in Fig. I A. II that there is a strong correlation between 
the "true" entropy ratio and the "pseudo" entropy ratios. We fit- 
ted the points in Fig. lA.ll with a linear model taking into account 
the errors on both coordinates (with the IDL Astrolib procedure 
FITEXY.PRC0) and we found 

— — = -0.051 +0.004 + (0.792 + 0.011)* (A.l) 

S OUT SQUT 



5 http://idlastro.gsfc.nasa.gov/ftp/pro/math/fitexy.pro 
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Table 1 : Observations analyzed for the clusters in the sample. In case of multiple observations, the one indicated with * (usually the 
longest) is the one used for spectral analysis. The exposure time in the third column is the sum of the exposure times of the three 
instruments after SP cleaning. In Col. 4 we list the symbols and colours used to mark each cluster in the figures of the present paper. 



Cluster 


Obs ID 


Exposure time (ks) 


Symbol 


A4038 


0204460101 


78.1 


Black open triangle 


A2199 


0008030201* 


42.8 


Black open square 




0008030301 


11.2 






0008030601 


11.3 




A496 


0135120201 


47.2 


Black open circle 


2A0335+096 


0147800201 


230.8 


Black filled triangle 


A2052 


0109920101 


85.1 


Black filled square 


A576 


0205070401* 


42.7 


Black filled circle 




0205070301 


28.1 




A3571 


0086950201 


43.5 


Blue open triangle 


A119 


0012440101 


54.1 


Blue open square 


MKW3s 


0109930101 


99.2 


Blue open circle 


A 1644 


0010420201 


42.0 


Blue filled triangle 


A4059 


0109950101* 


32.4 


Blue filled square 




0109950201 


64.9 




A3558 


0107260101 


126.4 


Blue filled circle 


A3562 


0105261301* 


116.6 


Red open triangle 




0105261501 


45.0 






0105261601 


52.5 






0105261701 


57.7 






0105261801 


23.6 




Triangulum Australis 


0093620101* 


27.3 


Red open square 




0093620201 


8.3 






0093620301 


31.8 




Hydra A 


0109980301 


52.2 


Red open circle 


A754 


0136740101* 


40.6 


Light blue open square 




0136740201 


14.8 






0112950401 


31.8 






0112950301 


34.8 




A3266 


0105260901* 


66.5 


Light blue open circle 




0105262001 


13.7 






0105262101 


16.2 






0105261101 


33.1 






0105262201 


9.1 






0105261001 


5.8 






0105260701 


56.4 






0105260801 


56.7 






0105262501 


17.3 




A85 


0065140101 


34.6 


Red filled triangle 


A3532 


0030140301 


27.7 


Green filled star 


A3667 


0206850101* 


162.9 


Light blue open triangle 




0105260101 


15.4 






0105260601 


62.5 






0105260401 


45.6 






0105260301 


47.2 






0105260501 


38.6 






0105260201 


31.8 




A2319 


0302150101 3 


27.6 b 


Red filled square 




03021 5020 l a 


28. 2 b 




A2256 


0141380201* 


33.1 


Light blue filled circle 




0112950601 


25.3 






0112951501 


28.1 






0112951601 


30.9 






0141380101 


24.4 




A3158 


0300210201* 


33.2 b 


Red filled circle 




0300211301 


14.0 b 




A 1795 


0097820101 


97.3 


Green open triangle 


A399 


0112260101 


27.9 


Green open square 
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Cluster 


Obs ID 


Exposure time (ks) 


Symbol 


A401 


0112260301 


34.7 


Green filled triangle 


A3112 


0105660101 


64.6 


Green filled square 


A2029 


0111270201 


30.8 


Green filled circle 


A2255 


0112260801 


25.1 


Pink open triangle 


A1650 


0093200101 


75.0 


Pink open square 


A2597 


0147330101 


144.3 


Pink open circle 


A478 


0109880101 


136.7 


Pink filled triangle 


PK.S0745-191 


0105870101 


34.3 


Pink filled square 


A2204 


0112230301* 


51.2 


Pink filled circle 




0306490101 


17.6 b 






0306490201 


30.9 b 






0306490301 


22.7 b 






0306490401 


29.7 b 




A1689 


0093030101 


106.7 


Pink filled star 



Notes: a We used both observations of A2319 to increase the statistics. b Since these observations have been performed after 
degradation of CCD 6 in MOS 1 in 2005, we did not consider the MOS 1 detector and the exposure time given in the table is just 
MOS2 + pn. 



